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In This Issue... 

Robert Hill’s Dithering As A Shading Technique 
reveals in detail the techniques used to generate 
the images and painterly backgrounds seen in the 
four slides in set #26, and in slide A2422 from set 
#27, distributed with this issue. 


Michael Freeman’s Mandelbrot And Julia Sets Of 
Elliptic Functions discusses the fractals that result 
from iteration over the Elliptic functions: a refresh- 
ing change from the quadratic, polynomial and 
rational functions most often used for iteration. 


In How To Spin M-sets On A Monitor Screen In Real 
Time, John MacManus discusses a program for 
rotating 3-D representations of Mandelbrot set 
Views. 


Dithering As A Shading Technique 
For Fractal Images 
Robert Hill 


Although there has been little mention of dithering 
in the literature on the display of Mandelbrot and 
Julia sets, this technique has much potential for 
producing beautifully shaded, often 3-dimensional 
appearing versions of these objects. This term, 
dithering, originally denoting a state of trembling 
or agitation, has come to mean in the art world the 
introduction of noise or randomness into the color- 
ing of various elements of a picture. This random- 
ness is a means of obtaining greater color fidelity 
at the expense of spatial fidelity or image sharp- 
ness. In representational art this is desirable to 
show the gradations in brightness of unevenly lit 
or shaded surfaces, for example. In the case of 
abstract computed images such as Mandelbrot sets, 
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where the color is arbitrary, dithering can create a 
similar effect of shading, thereby giving a physical, 
3-dimensional appearance to the mathematical 
object. In so doing, dithering replaces the discrete 
color steps defined by the Mandelbrot algorithm 
with a smooth appearing, sloping surface. 

There are many ways of dithering a picture of a 
Mandelbrot or Julia set. For any of them, however, 
the effect obtained is strongly dependent on the 
contrast in color and brightness of the adjacent col- 
ors involved in the dithering. Where the differ- 
ences in value and hue are small, the result will be 
a very smooth appearing surface. If the differences 
are larger, the surface will take on a grainy or 
speckled character. Slide A2421 of slide set S25 
affords a good example of this. The same dithering 
algorithm was applied to the entire picture, but the 
peripheral “background” gives an effect quite dif- 
ferent from that of the central figure. 

The very high contrast inherent in black and white 
pictures makes it impossible to demonstrate the 
more subtle effects seen with color. For this reason 
the reader is referred to the slide supplement for 
some examples. This importance of contrast also 
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makes it evident that in selecting a palette for a 

picture, the effect desired from the dithering must 

be taken into consideration. 

The simplest way to introduce dithering into a pic- 

ture is merely to add a random term to the expres- 

sion defining the color to be assigned to a pixel, as 
color = colorNumber + rnd*k 


For most purposes, k should be selected so that 
color often, but not always, assumes the next 
higher color number, or perhaps occasionally 
jumps to a still higher color. The greater the value 
of k, however, the greater the effect of obscuring 
the sharpness of the Mandelbrot image, so that val- 
ues greater than about 2 are not satisfactory. 


The Mandelbrot and Julia sets do have an advan- 
tage: outlying areas of low dwell are the regions 
where dithering is most effective. In the convo- 
luted fractal regions where contour lines are very 
close together, the variability inherent in the set 
overwhelms the dithering, so that the sharp edges 
are still retained where sharpness is most wanted. 


Conversely, in the outlying regions where contours 
are far apart, the contours tend to disrupt the dith- 
ering, and it is often desirable to obscure the con- 
tour lines as much as possible. One way to do this 
is to randomize the escape function, e.g. 


x+y <4+rnd. k. 


This has the effect of spreading the contour into a 
band of width k. In combination with the dithering 
routine, this technique can completely obscure the 
outlying contour lines while having little effect on 
the central figure. 


The Slides (S26, continued) 

Bob Hill’s four slides were sent out with Amy #26 
without technical data. Since then Bob has sent 
along the following information. The names of the 
slides are mine, not his. 


See also Bob’s article, Dithering as a shading tech- 
nique for fractal images, elsewhere in this issue. 


A2412: Fierce moth. A detail of the Julia set loca- 


ted at 0.2525 + 01, for the function Z + c. The tex- 
ture arises from dithering using a “checkerboard” 
formula of the form: 
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if px mod 2 + py mod 2 = 1 
then Color = D - 10*rnd 
else Color = D + 10*rnd 
where px, py = pixel coords, D = current dwell. 


was 400. Note that 
the dithering was applied only when D< 235. 


The number of iterations, D 


Max? 


A2417: Oils on rough canvas. A portion of the 
Julia set at 0.085 + 1.0077 of the function 

SS Së 
(z— zr + Sé - c, which is the power series for 
sin z truncated to three terms. The dithering for- 
mula was similar to that above, but was used in 
combination with a randomized escape limit func- 


tion (ELF), namely Ai +y <St+rnd-6. 


Note the texture resembling chalk on canvas; it is 
more conspicuous here, probably due to greater 
contrast between adjacent colors. 


A2413: Prickly pear. The central region of the Julia 
set at 1+ 0.17 for the same function as A2417, 
except that the two terms in the x, y expression 
have different divisors: 60 rather than 24. Thus the 
function differs somewhat from sinz. The dither- 
ing formula is also like that of A2417. A distinctive 
feature here is the ELF, with limit 16: it is the iter- 
ated function itself. This gives rise to the asymmet- 
ric swirls on either side of the central spiral. 


A2421: Spider's web. This is a detail (relative mag- 
nification 3.3) of Bob’s A2422. See below for a 
description of A2422. 


The Slides (S27) 


A903 (James Loyless): San Marco. The insides of 
Mandelbrot’s Dragon of San Marco are mapped 
using the minimum value of z as the inside color 
/ alternation criterion. Center =0 x 0.93, aspect 
ratio = 1.444, dwell limit = 250, Znin multiplier = 
250. 


A908 (James Loyless): Z-eye. An eye from Sea 
Horse land. This one is my wife’s favorite. 
Center = —0.74355+0.131995 x 58,800, aspect 
ratio = 1.41, dwell limit = 1000. 


A2422 (Bob Hill): Dragon. A Julia set of sinz 
(truncated) located at 1.13 + 0.4957. The dithering 
is similar to A2413, but with an added 0.015 - py 
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term (where py is the pixel coordinate), giving rise 
to the color gradation from top to bottom. The ELF 
is randomized as well. This calculation differs, 
however, in that the function z +c was used for 
the first two iterations. This slide is the “big pic- 
ture” of which A2421 is a detail. 


A4708 (Gérardin): PHILOSOPHER’S EMERALD. 
Julia; c = -1.476 (FRACTINT + spherical 
projection). 


Mandelbrot And Julia Sets Of Elliptic 
Functions 
Michael Freeman 


After seeing the beauty and complexity of Julia and 


Mandelbrot sets for the function ze 2° + c it is 
natural to wonder what we would see after iterat- 
ing other functions in a similar manner. Readers of 
Amygdala will be familiar with the results for other 
elementary functions, e.g. higher order polynomi- 
als, exponentials, and the circular functions sin z 
and cos z. Recently I have been exploring Julia 
sets of elliptic functions. Although these are less 
familiar than the functions mentioned above, they 
have some very nice properties and are not too dif- 
ficult to calculate. 


The theory of elliptic functions, developed in the 
19th century by Gauss, Jacobi, Weierstrass and 
other great mathematicians, is a wonderful accom- 
plishment of mathematical analysis [Abr64], 
[Whi27]. For the reader familiar with complex num- 
bers and elementary calculus I will attempt here to 
scratch the surface. (I apologize for any offence to 
serious mathematicians.) 


Examples of periodic functions are familiar to 
most of us: the circular functions sin and cos with 
the period 2r (radians). Specifically, this means 
that sin z= sin(z+27) for any complex number 
z. Less familiar at first glance are the hyperbolic 
functions sinh and cosh. These are defined in 

z =z 
-e 


2 H 


terms of exponentials: sinh z = 
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cosh z= . They are also periodic, but their 


period is along the imaginary axis, Amt Because 
the properties of the exponential function, we 
have sinh (ix) = —isinx, and cosh (ix) = cosx. 
The circular and hyperbolic functions are both 
periodic, each with a single period. 


Elliptic functions are doubly periodic. That is, they 
have two distinct periods in different directions in 
the complex plane. In fact, the general definition 
of an elliptic function is: any function of a complex 
variable which is doubly periodic and which has 
poles, but no essential singularities, in the finite 
plane. (The terms “pole” and “essential singularity” 
refer to points where the value of the function 
blows up or goes to infinity. The distinction is that 
the limit of the function as we approach a pole 
behaves like an inverse power of z. Near an essen- 
tial singularity, the behaviour is more compli- 
cated.) Another definition, apparently completely 
different, is in terms of the inverse of certain indef- 
inite integrals. One way to define the circular and 
hyperbolic functions is as the inverses of integrals 
of the square root of quadratic functions, e.g. 


d 
u Í a Ọ = sinu 
= —— s = H . 
eR 
This is generalized to elliptic functions by consid- 


ering integrals of the square root of cubic or quar- 
tic functions. 


I've investigated a type of elliptic function called 
Jacobian. These functions arc useful; they repre- 
sent the arc length of an ellipse and the displace- 
ment of a simple pendulum, among other things. 
The three commonest functions sn, cn, dn are 


d 
dt 

defined as follows. Given u we 

I — main t 

define sn(u) = sing, cn(u) = cos@, and 


dn(u) = J1— main 


The parameter m ranges from 0 to 1. To show the 
dependence on m we write sn(u|m) , cn(u|m) , 
dn(u|m). When m = 0, the elliptic functions 

reduce to circular functions. When m = 1, the ellip- 
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tic functions reduce to hyperbolic functions: 


sn(u|O) = sinu sn(u|0) = tanhu 


cn(u|0) = cosu cn(u|0) = sechu 


dn(w|0) = 1 dn(u|1(0)) = sech 


When m is between 0 and 1, these functions have 
periods along both real and imaginary axes. The 
periods depend on m: as m goes to 0, the real 
period goes to 2m and the imaginary period goes to 


ico, AS m goes to 1, the real period becomes œ and 
the imaginary period becomes 211. The periods 
can be found in terms of the complete elliptic inte- 
gral, wo), when @ = 2/2. This quantity is 
denoted K and depends on m. Also, K(1—m) is 
denoted K’. The real period is 4K and the imagi- 
nary period is 21K’, so 

sn(z) = sn(z+4XK) = sn(z+8K) =... 

= sn(zt+ 21K’) = sn(z+4ik’) = ... 


In addition, there are poles at z = iK’ and at 
2K+ ik’. 


These functions have many interesting properties, 
similar to the circular functions. For example, 


sn(—u) = —sn(u), can(—u) = cn(u) 
d e , 
Wyn = cn(u)- dn(u) 


(lm = sn(u4| 1- m) 
sn (iu = STEE EE 
sn(u + v) = 
sn(u)-cn(v)-dn(v) + sn(v)- cp. dn(u) 
1—m- en A0 . en Ta 
The first four of these properties follow in a 
straightforward way from the definition. The fifth is 


less obvious; see Whittaker and Watson [Whi27] for 
the proof. 


Calculating the elliptic functions 


I used the method of the arithmetic-geometric 
mean (agm — see appendix) which involves cal- 
culating trig, inverse trig, and exponential func- 


tions. Naturally, this takes more time than 2+ C; 
for speed I programmed my 83D87 math processor 
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in assembler which led to acceptable calculation 
times (25 MHz) using the Mariani algorithm. 


The actual functions iterated were: c-sn(z|m), 

c: cn(z| m), and c- dn(z| m), where c is a com- 
plex parameter. Since these functions are doubly 
periodic in z, the Julia picture is doubly periodic in 
zas well — good for wallpaper designs. The 
appearance is quite sensitive to the choice of 
escape radius (and of course, c and m). I made 
sure not to check for escape before the first itera- 
tion. 


Figure 1 shows the Julia set of 2.9 - cn(z| 0.2). The 


double periodicity can easily be seen, as well as 
the location of the poles. A maximum dwell of 100 


was used. The escape radius is ./20, as it is for the 
other figures as well. 


sé A er b CH WR KÉIS th P ii Ni D 


EC, E EE, 
a > 
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Figure 2. Mandelbrot Set of c- sn(z| 0.2) 


For Mandelbrot pictures, c is the screen point; but 
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one change is made from the usual procedure: z 
cannot be set to 0 for the sn functions, because 
sn(0) = 0. Instead, zis initialized to c. This works 


for Ze also, with one extra iteration. Figure 2 is 
the Mandelbrot set for c- sn(2 0.2) between —13.3 
and +13.3 on the real axis, and between —10 and 
10 on the imaginary axis. As that figure shows, 
there is a Mandelbrot-like entity in the centre, and 
the double periodicity is not exact. 


Figure 4. Central Feature of Figure 2, 
with Julia Sets 


A relationship between points on the Mandelbrot 
picture and the Julia sets exists. Figure 4 shows a 
blowup of the central feature of Figure 2 with four 
Julia sets (these are calculated between -2 and +2 
on the real axis, -1.5 and +1.5 on the imaginary 
axis. 


Figure 3. Detail of Figure 2 


For the elliptic functions, neither the Mandelbrot 
set nor the Julia sets from interior points are con- 


The dwell contours are very different from the 
usual sets. The window appearing in Figure 2 is 
blown up in Figure 3. Note the midget Mandelbrot. 
(For figures 2, 3, 4 the maximum dwell used was 
50.). 


nected. Nonetheless, it seems that there might be a 
property that would play the role here that con- 


nectedness plays for 2+. 


One could investigate these functions more 


deeply; so far I have just begun. 


Appendix: calculating sn, cn, dn 
When the parameter m is less than 0.05, we use: 
sn(u| m) = sin ( u) — 0.25m(u-— sin ( u) cos (u) ) cos (u) 
cn(u| m) = cos (u) + 0.25m(u-— sin (u) cos (u) ) sin (u) 
dn (u| m) = 1 — 0.5 msin? u 
When the parameter m is greater than 0.95, we use: 
sn(u| m) = tanh (u) + 0.25 (1 — m) (sinh (u) cosh (u) — u) sech? ( u) 
cn(u| m) = sech (u) — 0.25 (1 — m) (sinh (u) cosh (u) — u) tanh (u) sech (u) 
dn (u| m) = sech (u) +0.25 (1 — m) (sinh (u) cosh (u) + u) tanh (u) sech (u) 


b = 41- mM, CG = Jm, and iterate: 


Otherwise, we use the agm method: 


Given m, we create a table of triples (a, D.C BD sg 
for n = 0, 1, 2, ..., N. We start with a, = 1, an = 2 
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So 
3 
Il 
x 
i 
S~ 
a 


until Cy is less than 0.001. ON will not exceed 12.) 


We then create a table of om, starting with 


Qy = ae Gy: u, and computing the rest from the 


C 
relation sin(29,_,-@,) = z Sn (o 1. 


Finally, we have 
sn(u| m) = sin (Q,) 
cn(u| m) = cos (,) 
COS (Q,) 


dn (u| m) = cos (@,— 9) 


Now, we require the relationship between the 
functions of a complex z and those of its real and 


imaginary parts (x, y). These are: 
S:a,+icd-Ss,-¢, 


sn(x+ iy) = Se 
CG (SO: D: d, 
cn(x+ iy) = aan ; aie 
d-c,:d,—imsc:s, 
dn(x+ iy) = eos; aia 


where D = Ee mS sS, s= sn(x m), 
C= cn(x| m), d = dn(x| m), s, = sn} 1- m), 
c& = cn(y| 1- m), and d, = dn(y|1 -— m). 
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Corny M 


The world’s most mathematical hoaxers struck in 

the cornfields of Cambridge last week. Rather than 
a mere corn circle, they created a corn Mandelbrot 
set — the most complicated object in mathematics. 


Corn circle of the chaotic kind 
The possibility that intelligent aliens have mastered 
fractal geometry just a couple of years after 
humans managed it seems unlikely. Instead suspi- 
cion has fallen on the University of Cambridge 
mathematics department. “We know nothing about 
it,” claimed Peter Landeshoff. Benoit Mandelbrot, 
the IBM researcher who discovered the set, also 
denies any connection. He says he was in North 
America all last week. 


The image went unrecognised by the Cambridge 
Weekly News which published this picture last 
Thursday. On Friday the paper printed a list of 
readers who pointed out that it probably owed 
more to Unemployed Fractal Observers than the 
other kind of UFOs. Unlike traditional shapes, frac- 
tals like the Mandelbrot set cannot be drawn with a 
pencil, gusts of wind or a scythe. They have to be 
built up dot by dot. However, mathematicians have 
recently proved that it is impossible to draw the 
Mandelbrot set accurately, even dot by dot. A 
really intelligent alien would have known that. 


Mandelbrot munchies 

In regard to the Mandelbrot corn-circle, you cor- 
rectly state that it is impossible to draw the Man- 
delbrot set accurately dot by dot, as the images 
would be constructed in a cornfield. This would 
suggest we look ever closer at the implications of 
cornfields. 
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Within the swirls and seahorses of the set, once 
iterated, small but perfect Mandelbrot sets present 
themselves. Where should we look for these sub- 
sets? 

I am pleased to be able to enclose the first piece of 
sold evidence towards a proof of the hypothesis 
that a higher intelligence is attempting to establish 
fractal communication with us. You will note that 
the ridge near the point of the main set visible on 
the cornflake consists of a series of peaks that may 
well be worth further examination under a scan- 
ning electron microscope. 


Mandelbrot pop 

Regarding the Mandelbrot cornflake, it follows that 
if a Mandelbrot cornflake is the self-same at all 
scales, then the Universe originated not with a big 
bang but a snap, crackle and pop. 


Reprinted from New Scientist (24 August 1991) p.16 ff. 


How To Spin ™-sets On A Monitor 
Screen In Real Time 
John MacManus 


A while back I came across an ad for a program 


called AcroSpin! which promised to display 3D 
objects on a PC screen in any video mode. In addi- 
tion, the program promised to permit the normal 
manipulations of scaling, translating, and rotating 
in real time. Believe it! 


The input data file for AcroSpin is an ASCII list of 
XYZ coordinates, and simple line, point, and color 
statements. The input file can be written in any 
vanilla text editor. My first input file was written by 
hand in my word processor with instructions to 
join some 350 atoms of a protein molecule in 3- 
space. It was worth it. The result was a protein 
backbone rotating smoothly on my home AT. Now 
I use semiautomatic ways to translate protein coor- 
dinates from a protein data base into AcroSpin for- 
mat. 


However, my interest in graphics is also a hobby, 
and so blossomed an interest in adventures in the 
complex plane. I looked at all the pretty pictures in 


1. From Acrobits, PO Box 5563, Redwood City, CA 
94063-0563; $30.00. 
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the shiny books and drooled. I worked my way 
through the code in Becker and D6rfler’s book 


Dynamical Systems and Fractals’ and grew my 
own fractals at home. I was particularly interested 
in getting as much information onto the screen as 
possible, and so was particularly attracted to the 
pseudo-3D plots produced by incrementally offset- 
ting the X-dimension to the right. I assigned the 
iteration number to them third dimension, and 
plotted it in the ¥ dimension on the screen. This is 
done very cleverly by Becker and Dörfler from 
front to back with hidden line removal. An alter- 
nate approach by Ken Phillip with similar results 
recently came to my notice in a back issue of AMY 
(#13, October 1988). 


But why not real 3D images? Why not output M- 
set coordinates to AcroSpin using the iteration 
number as the Z dimension? Then one could rotate 
the fractal to one’s heart’s content. AcroSpin has a 
limit of about 10,000 points according to the man- 
ual. Since an EGA screen has 224,000 points, a full 
screen image is out of the question, but how about 
limiting oneself to a quarter EGA screen (56,000 
points), and plotting every other pixel in both 
directions? That would give a usable file of 14,000 
points. However, the data file is over 100 Kbytes, 
and the images are very sluggish to move because 
of the large block of points inside the M-set. 


How do we lighten up the image? This can be 
accomplished by plotting only a few points on the 
base plane to indicate its horizontal dimension. In 
addition, once inside the M-set, instead of drawing 
each point, lines may be drawn between Rz and 
Xma, for any given Y. This reduces data files to 40 
- 60KB, which will rotate beautifully. AcroSpin 
sucks the entire data file into memory, and auto- 
matically scales the final image to be nearly full 
screen. The image can be scaled down to fill in the 
holes used to decrease data size. It can be viewed 
in whatever video mode (monochrome or color) is 
available to the user, including CGA, EGA, or VGA. 


Another interesting feature of AcroSpin is layering 
of the image. In this case, ranges of iteration num- 
ber can be assigned not color only, but also differ- 
ent layers of the image. This permits dissection of 


2. Cambridge University Press, 1989) 
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the image’s iteration profile, and permits analysis if 
only monochrome is available. 


The program to produce output in AcroSpin format 
was written in MicroSoft’s QuickBASIC. The M-set 
algorithm is the usual one familiar to Amy readers, 
and consists of a mapping technique along the X- 
axis with incremental Y including the common test 
for inside or outside the set. If the point is outside 
the M-set, then the Z dimension is assigned the 
iteration number. The X, Y and Z are written to 
disk with the AcroSpin point notation. If the point 
is inside the set, then that point is flagged in an 
array sized to the X dimension. Before increment- 
ing the Y dimension, the array is reexamined for 
beginnings and ends of the set boundary. The 
boundary X points (for later line drawing), Yand Z 
(the latter being the inputted maximal iteration) 
values are plotted to disk. Alterations in coordi- 
nates, iteration number, color, or layer can be 
made at the user’s discretion. 


The QuickBASIC listing to produce output of M-set 
coordinates suitable for display by AcroSpin can be 
obtained for postage cost from: John D MacManus 
/ 38 Fentiman Ave. / Ottawa, Ontario, K1S 0T6 / 
CANADA. 


Effects Of The Shape Of The Escape 
Region On ™-set Images 


(Hey Charley—where’s the T-shirt?) 
Ken Philip 


In AMYGDALA #9, Charles Fitch presented an 
image of a region in west Sea Horse Valley that 
showed a very odd structure consisting of con- 
nected lines of black and white blobs that ran per- 
pendicular to the normal contours seen in 
Mandelbrot Set images. There was also a color 
slide showing the same region, which showed a 
progression of color shades along the rows of 
blobs — suggesting that the connected black (or 
white) blobs did not have the same dwell. He then 
offered to provide a T-shirt to the first person who 
could discover the boundary shape for the escape 
region used to create these images. 


I had not made a determined effort to solve this 
puzzle until recently, when a question about the 
influence of the shape of the escape region on the 
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Figure 5. Fitch’s puzzle picture 
Alternate dwells black and white 


Escape region 2 + 10.,/1—|cos28| 


resulting image turned up on the FracL-L confer- 
ence on BITnet. That brought Fitch’s puzzle to my 
attention again, and I began modifying a True 
BASIC Mandelbrot program in an attempt to repro- 
duce his figure. I started out by looking into wrap- 
ping trig functions around the origin, letting the 
escape ‘radius’ be a sine or cosine function of the 
angle. I got some interesting-looking images, but 
nothing like Fitch’s until I made a mistake and 
ended up with a large coefficient on a cosine func- 
tion. Figure 5 shows the result, and indicates that I 
was on the track of the solution, if not very near it. 
Although the shapes of the black blobs were noth- 
ing like those in Fitch’s image, the blobs were 
arranged along rows perpendicular to the dwell 
bands. I realized that the important thing here was 
that the value of root quite large at certain angles 
— and looked for other functions that had the 
same property. 


Figure 6 (next page) shows the image for a func- 
tion utilizing the tangent of the angle. The differ- 
ences between this and Fitch’s image are rather 
small. 


I had been graphing the trig function in polar coor- 
dinates to see what the escape region looked like, 
and it became apparent that the best match to 
Fitch’s image was produced by an escape region 
that was approaching two vertical lines symmetri- 
cally placed around the origin, as indicated in Fig- 
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Figure 8. Fitch’s puzzle picture; True BASIC Man- 


Figure 6. Fitch’s puzzle picture delbrot program; Alternate dwells black and white; 
Escape region —2 < R&c) <2 


Alternate dwells black and white 
Escape region 3 + Jang 
Figure: 7. 


The dotted circles about the origin in Figure 7 are 
the ‘normal’ escape radius r= 2. At this point the 
light bulb finally turned on. Since I suspected that 
Fitch had used something simple, I tried testing for 
the real part of Z going less than -2 or greater than 
+2, while ignoring the imaginary part of Z (see the 
right hand image in Figure 7). 


The result was Figure 8, which is a precise match 
of Fitch’s figure in AMYGDALA #9. 


Figure 9 shows a MandelZot image of the same 
area for comparison. 


So the puzzle was solved, if by a rather round- Figure 9. Fitch’s puzzle picture 
about route... MandelZot image, dwell-band width = 100 


Escape region for Figure 5 Escape region for Figure 6 Escape region for Figure 8 


I should add that, as inspection of Figure 5, Figure 
6 and Figure 7 will show, one can produce images 
of the M-set when the escape radius is set to less 
than 2, and the images can be reasonably well- 
behaved provided the region of interest is itself 
well inside of the escape radius. 


More On Fractal WitchCraft and 
Synchronous Orbits 

Steve Stoft asked for space to reply to Leonard 
Herzmark’s letter which appeared in Amy #26 as A 
Review Of Fractal WitchCraft. ’m reluctantly 
acceding to that request — “reluctantly” because I 
hesitate to perpetuate what may be moving 
towards an unseemly spat. “Acceding” because I 
guess I feel that Steve should have the right to 
reply. 

Also, I apologize for leaving the title of Leonard’s 
letter out of the Table of Contents in Amy #26. 


Leonard Herzmark was twice asked to review FW 
(Fractal WitchCraft). On the first occasion he com- 
plained that none of our list of commonly used 
commands worked. After a phone call and several 
days of scratching my head, I finally deduced that 
Leonard had ignored the only instruction at the top 
of the list: “Use these from the Main Menu.” He 
had just used them from wherever he happened 
be: 


So we added a “/” command to FW that returns 
you from anywhere to the main menu, and began 
every command sequence in the list with a slash. 
This time Leonard had no problem. 


Well, none till he tried entering coordinates: “of 
course not being able to enter the Ymin, I did not 
get the same graphics.” Fortunately he sent me his 
four “NOT” maps with the wrong images. Looking 
at his entered coordinates I found what I sus- 
pected; Leonard had entered the Ymin value into 
the Ymax input slot. He had then wanted to enter 
Ymax into what he correctly saw to be the Ymin 
slot, but as FW tells you, it just calculates this value 
form the other three. Leonard had no problem with 
FW, he was just unable to figure out which of the 
two positive Y values given in B-of-F was the big- 
ger number! Entering Ymax in the Ymax slot gives 
the right graphic every time. 
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If you have trouble telling which of two numbers 
is bigger you may still want to purchase FW for the 
trippy slide show. @ 


Bert Tyler, one of the co-authors of Fractint, I 
believe, sent Steve the following assessment of FW: 


I'll confirm that FW’s algorithm appears to be sig- 
nificantly faster than Fractint’s when used with 
images that require high maxiter values. It would 
appear that you folks have found another way to 
beat brute force with intelligent algorithms. Con- 
gratulations! 


Book Review 
Review by lan Entwistle 


Fractal Geometry — Mathematical Foundations 
and Applications, Kenneth Falconer. Wiley 
(Chichester & New York, 1990), hard cover, 310 
pp. $19.95 


The contents page of the book lists eighteen chap- 
ters. Chapter headings range from Hausdorff mea- 
sure and dimension, Dynamical systems, Projec- 
tion of fractals, Iteration of complex functions — 
Julia sets, to Random fractals and Multifractal 
measures. 


Most topics are treated by a section on background 
followed by theorems or propositions. Many exam- 
ples follow. Each chapter is ended with a list of 
CXErcises: 


Although the consequences of many of the proofs 
may be of interest to the general fractal enthusiast, 
this is a textbook for mathematicians or mathemat- 
ics students. 


Many of the diagrams, which include Koch curves, 
coastlines, Julia sets, and various attractors, will be 
familiar to most casual readers. The accompanying 
maths text will not be so easily followed. 


The substantial reference list (alphabetical by 
author) is up to date. 


This book will convince most mathematicians 
unfamiliar with fractal geometry that fractals are 
not just the pretty pictures being generated by 
computers to please the media and pattern genera- 
tors. They are seen to embody complex and fasci- 
nating mathematics. This book should be all a 
mathematician needs for conviction. 
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In the generally more readable introductory chap- 
ter, the author adopts a more rigid view about frac- 
tal than has developed in other recent texts. His 
detailed mathematical treatise therefore carefully 
avoids the wider application of fractals declared by 
B. Mandelbrot. 


For Amygdala readers who are fractal generators 
the book cover illustration will be of interest. A 
nice coloured rendering of the Julia set for 


2 + 0.233 + 0.53787 certainly attracts the eye. 


This book is not recommended unless you are a 
trained mathematician. 


Author 


— Mike Freeman 


Born and lived almost all my life in Vancouver, 
B.C. Degrees in Physics from UBC and Cal Tech. 
My primary occupation is teaching physics at 
Capilano College (2 yr.), but for the last 10+ years 
I've been doing programming and systems analysis 
for the college. My attraction to the Mandelbrot Set 
is twofold: the fascination with the set itself, and 
the challenge of writing better and faster programs 
to display it. 
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Books 
— from the Springer Newsletter, Computer Science 
#1/92 


Fractals for the Classroom, H.-O. Peitgen, H. Jür- 
gens, D. Saupe, et. al. 

Four new books on fractals and chaos. The 
authors, with help from leaders of the teaching 
community, have created books which are setting 
standards for science and teaching. The well writ- 
ten, clearly illustrated texts explore both the history 
and the unlimited potentials of fractals. Activities 
books encourage participation with exercises that 
are designed to develop a basic mathematical 
understanding and appreciation of fractals. 


The first two volumes consist of independent chap- 
ters covering the major concepts related to the sci- 
ence and mathematics of fractals. Written at the 
mathematical level of an advanced secondary or 
beginning undergraduate student, Fractals for the 
Classroom, Parts I & IJ include many fascinating 
insights for the instructor, and integrate illustrations 
from a wide variety of applications with an enjoy- 
able text to bring the concepts alive. These texts 
will have quite an impact upon professors, stu- 
dents, and the mathematics education of the gen- 
eral public. 


The second two volumes discuss strategic activities 
which are designed to develop, through a hands- 
on approach, a basic mathematical understanding 
and appreciation of fractals. The exercises have 
been developed from a sound instructional base, 
stressing the connections to the contemporary cur- 
riculum as recommended in the National Council 
of Teachers of Mathematics’ Curriculum and Evalu- 
ation Standards for School Mathematics. The slide 
package that accompanies the first volume 
includes some of the highest quality fractal images 
available anywhere. 


Introduction to Fractals and Chaos; 450 pp., 289 
illus., hardcover $29.00 (1992). 


Complex Systems and the Mandelbrot Set; 384 pp., 
141 illus., hardcover $29.00 (1992). 


Strategic Activities, Volume I; 128 pp., 143 illus., 9 
color slides/ Softcover $19.95 (1992). 


Strategic Activities, Volume II; 187 pp., 75 illus./ 
Softcover $19.95 (1992). 
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The Science of Fractal Images, HO. Peitgen and 
D. Saupe. 

“Anyone interested in the internal logic of the frac- 
tal world will find the book irresistible.” 

— Nature 


Offering a comprehensive survey of the fractals 
that are most interesting from a graphics perspec- 
tive, this “how-to” book demonstrates the way to 
create these striking images on a computer. Con- 
tributors include M.F. Barnsley, R.L Devaney, B.B. 
Mandelbrot, H.-O. Peitgen, D. Saupe, and R.F. 
Voss, with additional contributions by T. Fisher 
and M. McGuire. 


312 pp., 142 illus., 39 color plates/ Hardcover 
$39.95 (1988). 


Fractals and Chaos, T. Crilly, R.A. Earnshaw, and 
H. Jones. 

What are the relationships between the two 
domains: chaos and fractals? This volume brings 
together a number of distinctive contributions con- 
cerning this central issue. The book provides an 
understanding of fractals, chaos, and the interrela- 
tionship between the two domains, leading to a 
common basis for examining the growth, develop- 
ment, organization, and behavior of complex 
dynamic systems. 


277 pp., 146 illus./ Hardcover $39.50 (1991). 


The Algorithmic Beauty of Plants, P Prusinkiewicz 
and A. Lindenmayer. 

This book explores mathematical models of the 
developmental processes and structures of plants, 
and illustrates them by using state-of-the-art com- 
puter-generated images. The authors outline the 
modeling and rendering of plants which are suit- 
able for realistic image synthesis; the scientific 
potential of computer graphics in the visualization 
of biological structures and processes; the relation- 
ship between control mechanisms employed by 
living plants and the resulting complex develop- 
mental sequences and structures, and the relation- 
ship between developmental processes, self- 
similaritv, and fractals. 


228 pp., 150 illus., color plates/ Hardcover $39.95 
(1990), 


e The software and all the books are available 
from Amygdala (see Price List / Order Form). 
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Fading Away? 


Well, we’re down to 528 subscribers. Are people 
losing interest in fractals? Is it because of the long 
hiatus in publication over the past year? Has this 
newsletter become irrelevant or boring? Those 
questions could only be answered by taking a mar- 
ket survey, if at all, so Pll regard them as unan- 
swerable, except by each of me/you for him/ 
herself. 


I'll assume that many of you are interested in keep- 
ing Amygdala going. For that to happen I’m going 
to need your help to boost circulation. I believe 
there are several thousand potential subscribers 
out there; the problem is, how to reach them. I’m 
not interested in “selling” them, just in getting them 
to know about Amygdala. Here’s what you can do: 


1. Get fractal friends to subscribe. 


2. If you publish a book or periodical, think about 
mentioning Amygdala in it. 


3. If you are active on electronic forums or bulletin 
boards where fractals are of interest, please men- 
tion Amygdala there. 


4, If you read publications in which fractals are a 
topic, send a letter mentioning Amygdala. 


5. If you’re sharing your issues with someone else, 
get them to subscribe on their own. Saving money 
by piggybacking may be appropriate for big, rich 
publications like the Wall Street Journal, but it may 
just serve to kill the Amygdalan goose. 


Circulation 


As of 27 July 1992 Amygdala has 528 subscribers, 
167 of whom have the supplemental color slide 
subscription. 


Renewal 

For 33 of you subscribers out there this is the last 
issue of your current subscription. You can deter- 
mine if you are one of those by looking at the 
lower left of your mailing label: If it reads “A27” 
this is your last issue. if it reads “S27” this is your 
last slide set. 

I urge you to use the enclosed form to renew your 
subscription promptly to avoid missing anything. 
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TOTAL WEIGHT (Lines [2]+[4] from other side: pounds) i.is [7] 
For countries other than U.S. or Canada, use the following table to find your Rate Group: 


Rate Group Countries 


A Bahamas, Belize, Bermuda, Jamaica, Mexico. 
B Austria, Bolivia, Colombia, Luxembourg, Poland, Switzerland. 
C Denmark, Germany, Great Britain, Greece, Hungary, Ireland, Israel, Italy, Netherlands, 
Portugal, Spain, Taiwan, Turkey, Yugoslavia. 
D Algeria, Argentina, Australia, Belgium, China, Finland, India, Malaysia, New Zealand, 
Norway, Philippines, Sweden. 
E Brazil, France, Japan, USSR. 
U.S. Book Rates: For books ONLY to be mailed to an address in the U.S. ONLY, you can use 
the following special 4th class book rate (use weight from line [4]: 
1# Rate 2# Rate 3# Rate 4# Rate 5# Rate Add’l # 
Ath cl: FOS 1.48 1:91] 2.34 277 0.43 
POSTAGE & HANDLING, BOOKS ONLY (U.S. 4th Class Book Rate? [8] 


Air Parcel Post Rates (Priority Mail in U.S.) 


Use the following table to compute Postage & Handling (use weight from line [2] if you used the 
U.S. Book Rate for books; otherwise use weight from line [7]: 


1# Rate 2# Rate 3# Rate 4# Rate 5# Rate Add’! # 


LS: 2.90 2.90 = AIO 4.65 5.45 Consult postoffice 

Canada 5.00 5.00 6.40 7.80 9.20 1.40 

A 6.00 9.00 12.00 15.00 18.00 2.00 

B Jato 12.00 16.25 20:50 24.75 3.00 

$ Se 14.25 19220 24.25 29.25 4.00 

D 10.70 16.70 22.70 28.70 34.70 5.00 

E 12.30 19.30 26.30 ILAG 40.30 6.00 

POSTAGE & HANDLING et e e a el te A E [9] 
BOOKS / PORTFOLIO / VIDEOTAPE SUBTOTAL ([6] from other aide? [10] 


For shipments outside the U.S., compute insurance from the following table, based 
on the value from line [10], rounded up: 


Insurance Rates 


Value > $50 $100 $200 $300 $400 $500 

U.S. Insurance is included in prices. 

Canada 75 1.60 2.40 3.50 4.60 5.40 

elsewhere 1.60 2.40 3.50 4.60 5.40 6.20 

UNSURA e EE [11] 
NEWSLETTERS / SLIDES SUBTOTAL ([1] from other side? [12] 
TOPAL Si 19) 16) eo ATD EE [13] 
Payment enclosed (please check one): check J MasterCd hd VISA A] Other (specify) 

Credit Card #: Expiration Date: 

Name: 

Address: 


City, State, Zip: 

Signature (needed only for credit card): 

Telephone: 

Please pay in U.S. funds, with any of the following: A check payable through 


a U.S. bank, money order, postal money order, bank draft, cash (mailed at 


your own risk), credit card (VISA or MasterCard). July 27, 1992 


